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Abstract 

Minimal Dark Matter (MDM) stands as one of the simplest dark matter sce¬ 
narios. In MDM models, annihilation and co-annihilation processes among the 
members of the MDM multiplet are usually very efficient, pushing the dark matter 
mass above 0(10) TeV in order to reproduce the observed dark matter relic den¬ 
sity. Motivated by this little drawback, in this paper we consider an extension of the 
MDM scenario by three right-handed neutrinos. Two specific choices for the MDM 
multiplet are studied: a fermionic SU{2)l quintuplet and a scalar SU{2)l septu¬ 
ple!. The lightest right-handed neutrino, with tiny Yukawa couplings, never reaches 
thermal equilibrium in the early universe and is produced by freeze-in. This creates 
a link between dark matter and neutrino physics: dark matter can be non-thermally 
produced by the decay of the lightest right-handed neutrino after freeze-out, allow¬ 
ing to lower significantly the dark matter mass. We discuss the phenomenology of 
the non-thermally produced MDM and, taking into account significant Sommerfeld 
corrections, we find that the dark matter mass must have some specific values in 
order not to be in conflict with the current bounds from gamma-ray observations. 
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1 Introduction 


The identity of the Dark Matter (DM) that makes up about 25% of the energy content 
of the universe is one of the most important open problems in (astro-)particle physics. 
Lots of candidates have been proposed under the assumption that DM is made of parti¬ 
cles. The most popular options include Weakly Interacting Massive Particles (WIMPs), 
axions, gravitinos and asymmetric dark matter. In particular, WIMPs are theoretically 
well-motivated candidates since the present relic density of DM can be reproduced by 
thermal production in the early universe with an electroweak scale DM mass and an anni¬ 
hilation cross section (cm) 3 x 10“^® cm^/s, in the typical range for a particle with weak 
interactions. This intriguing coincidence, usually called the WIMP miracle, has triggered 
a massive exploration of WIMP DM scenarios, with detailed studies of their phenomeno¬ 
logical implications and dedicated experimental searches for DM in the form of WIMPs 
in colliders as well as in direct and indirect detection experiments. 

In scenarios with electroweak scale DM, a discrete symmetry is often imposed in order 
to stabilize the DM candidate. Although this symmetry is usually introduced by hand, 
many theoretical justihcations are known. For instance, this symmetry can be seen as 
a remnant after the spontaneous breaking of a larger symmetry group. Many proposals 
in this direction exist, based on global [1-3] or gauge symmetries [4-6], some of them 
linked to flavor symmetries [7,8]. A completely different approach is to consider that the 
origin of this symmetry is accidental. If a large multiplet of the SU{2)l gauge symmetry 
of the Standard Model (SM) is introduced, an accidental Z 2 symmetry may appear due 
to the restrictions imposed by gauge invariance and renormalizability. This is the so- 
called Minimal Dark Matter (MDM) approach [9], a popular scenario with some recent 
works [10-14]. 

There is, however, a generic drawback in MDM scenarios: the components of the large 
SU{2)l multiplets are generally required to be nearly degenerate. The origin of this mass 
degeneracy is easy to understand. First of all, in some MDM models this is actually 
predicted, since the mass splittings among components of the large SU{2)l multiplets 
only appear at the one-loop level [15,16]. When this is not the case, and large mass 
splittings can in principle be obtained, one must face two problems. First, the large 
SU (2 )L multiplets contribute to electroweak precision observables (EWPO) through the 
STU oblique parameters and, in order to be compatible with the current experimental 
measurements, one typically requires small mass splittings [17, 18]. And second, the 
mass splittings must be induced by Higgs-DM-DM couplings for scalar DM, which are 
required to be small in order to suppress the elastic scattering with nuclei via Higgs 
exchange and be compatible with direct detection constraints (see for example [19]). As 
a result of this degeneracy, all members of the multiplet will be in thermal equilibrium 
during freeze-out, co-annihilating very efficiently with the DM particles and strongly 
reducing the DM relic density. In order to reproduce the relic density measured by Planck, 
= 0.1186 ± 0.0020 at 68% Conhdence Level (CL) [20], this implies a heavy DM 


2 


particle. In fact, once Sommerfeld corrections are included [21], the DM mass is typically 
found to be about O{10) TeV. Although perfectly plausible, this is not very attractive 
from a phenomenological point of view. 

Another open problem that calls for physics beyond the SM is the existence of non-zero 
neutrino masses and mixings, as established by neutrino oscillation experiments. Many 
extensions of the SM have been proposed to address this issue. The energy scale for the 
new states responsible for neutrino mass generation can be either very high, potentially 
relating neutrino masses to unihcation, or low (TeV scale or below). In the latter case, 
the existence of new physics at low energies leads to interesting phenomenological per¬ 
spectives, within the reach of current collider and low-energy experiments. Furthermore, 
many neutrino mass models include DM candidates, although an interplay between these 
two fundamental issues is often missing. 

In this paper, we consider a very economical extension of the MDM scenario. In ad¬ 
dition to the multiplet containing the DM particle, three right-handed neutrino singlets 
are introduced. No additional symmetry for DM stabilization is required due to an acci¬ 
dental Z 2 symmetry resulting from the gauge invariant renormalizable interactions of the 
DM multiplet. The right-handed neutrinos play two roles in this model. First, neutrino 
masses are induced at tree-level with the standard Type-I Seesaw mechanism [22-26], and 
second, the out-of-equilibrium decay of the right-handed neutrinos leads to non-thermal 
production of DM, allowing one to compensate the strong effect of co-annihilations and 
lower the DM mass below the TeV scale.^ This setup will be illustrated with two spe- 
cihc examples: a model with a fermionic quintuplet and a model with a scalar septuple!, 
in both cases with vanishing hypercharge.^ The number of physical parameters in both 
models is limited and many experimental constraints exist. As a result of this, the setup 
is very predictive and will dehnitely be tested in future experiments. 

The rest of the paper is organized as follows: in Sec. 2 we introduce the setup, comment 
on the degeneracy among the members of the MDM multiplet and discuss the mechanism 
responsible for neutrino mass generation. The DM phenomenology of these models is 
explored in Sec. 3, where we present the main results of this paper. The most relevant 
constraints in our scenario are discussed in Sec. 4 and, hnally, we summarize and present 
the main conclusions of the paper in Sec. 5. 

2 The Models 

2.1 New particles and interactions 

We consider an extension of the SM by three right-handed neutrinos W (z = 1 — 3) 
and a SU{2)l multiplet y which is either a quintuplet fermion or a septuple! scalar, as 

^Non-thermal production of Wino DM (triplet) has been discussed in Ref. [27]. 

^Septuplets with hypercharge Y = 2 have also been considered in Refs. [28,29], although with a 
completely different motivation: the extended model keeps the p-parameter as 1 at tree-level. 
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iV* 

X 

Ff/(2)i 

1 

5(7) 

U(l)y 

0 

0 

spin 

1/2 

1/2 ( 0 ) 


Table 1: New particle content and charge assignment in the two models under consid¬ 
eration. Here i = 1,2, 3. The hrst model introduces a quintuplet fermion, whereas the 
second introduces a septuplet scalar, in both cases with Y = 0. 


summarized in Tab. 1. The quintuplet fermion and the septuplet scalar can be denoted 
by the vectors 



/ +X+++ \ 

+X++ \ 


+X++ 

+X+ 


+X+ 

-x° 

, X = i 

-x° 

+x" 


-X" 

+x“ y 


+x"" 


\ -X / 


The prefactor i and the sign for each component are taken in order to satisfy = x 
where x'^ denote charge conjugation of the held x- Further details about products of 
SU{2) multiplets are given in Appendix A. The kinetic and Yukawa Lagrangians of the 
new particles are given by 


C 


/-MDM 


+ - mjvJ Ni 


{yl^NiPLl^ + H.c.) , 


( 2 ) 


where is the kinetic term of the multiplet given by 


/-MDM 


-x^ {ip-m^)x 

i (D^x) (D^x) 


for the quintuplet fermion, 
for the septuplet scalar. 


( 3 ) 


The covariant derivative is dehned by = d^ + ig 2 W^T^ with the 811(2)^ gauge coupling 
g 2 , gauge held kF“ and generator T“ (a = 1,2,3). Note that Eq. (2) is written in the 
basis in which the right-handed neutrino Majorana mass term is diagonal. From the 
Lagrangian, one can see that the right-handed neutrinos dominantly interact with the 
longitudinal mode of gauge bosons and Z, while the MDM multiplet x interacts with 
the transverse mode. The scalar potential for the model with the quintuplet fermion is 
exactly the same as in the SM (V5 = Vsm), whereas that for the septuplet scalar model is 


F? — —+ 


T 




k=l 


+ 


Vx I 




( 4 ) 
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where 0 is the SM Higgs doublet with hypercharge = 1/2. In the following, the scalar 
coupling is assumed to be small, as in Ref. [9]. 

Several comments about the scalar potential in Eq. (4) are in order. First, we note that 
in some cases various singlet contractions of SU{2)l higher dimensional representations 
are possible. For example, four kinds of singlets can be constructed for the term. 
However only two of them are independent. This is denoted with a summation in Eq. (4). 
Furthermore, there is no septuplet cubic term since the singlet obtained by contracting 
three septuplets is completely anti-symmetric: 7®7(8)7 d1^. As a result, an accidental 
Z 2 symmetry appears, under which the septuplet y is odd and the rest of the fields even. 
The same accidental symmetry appears for the quintuplet fermion as well. Therefore, the 
lightest state contained in the quintuplet or septuplet will be stable and a DM candidate 
is automatically included in the model without any additional symmetry.^ One should 
note that for the septuplet scalar, the dimension 5 operator y^|0p/A, where A is the 
energy scale at which the operator is induced, leads to the decay of the DM candidate at 
one-loop level and is not forbidden by any symmetry. As it has been recently discussed 
in Ref. [30], A > 10^° GeV is necessary for ~ 10 TeV in order to obtain a DM lifetime 
longer than the age of the universe, tu ~ 10^® s. Moreover, pairs of gauge bosons 77, 'jZ, 
ZZ and W~^W~ are produced by the DM decays, and these are constrained by gamma- 
ray experiments. One finds that the DM lifetime should be roughly > 10^^ s to evade 
them. This corresponds to A > 10^^ GeV which is much larger than the Planck scale. In 
contrast, there is no possible 5 dimensional operator for the quintuplet fermion case [30]. 
The possible 6 dimensional operators are and where is 

the held strength tensor for the SU (2 )l gauge group. For such 6 dimensional operators, 
a DM lifetime long enough to satisfy the gamma-ray constraints can be achieved with 
A = 10^^ GeV. Thus, the quintuplet fermion DM candidate would be stable even if one 
includes non-renormalizable operators. 

Finally, lepton number is softly broken by the Majorana mass term of the right-handed 
neutrinos. As a consequence, neutrinos acquire Majorana masses through the canonical 
type-I seesaw mechanism, as we will see below. 

2.2 Mass degeneracy within the MDM multiplets 

After electroweak symmetry breaking, only the SM Higgs doublet gets a non-zero 
VEV, 0° = (0) -|- h/\/2. The Higgs boson mass is given by = A^(0)^ and all the 
components of the multiplet y have the same mass at tree-level. The mass is 
given by the bare mass in Eq. (3) for the quintuplet fermion and ~ /^x 
the septuplet scalar. However, a mass difference is induced at the one-loop level. Given 
two components of the MDM multiplets with electric charges Q and Q', the one-loop 

^For the septuplet scalar, we will assume that the vacuum of the theory is such that no other scalar 
held besides the usual Higgs doublet has a vacuum expectation value (VEV). This guarantees that the 
accidental Z 2 symmetry remains after electroweak symmetry breaking. 


5 



induced mass splitting is computed to be 


rriQ - mQ> = - Q'^) ^ 


Otw 


\mxJ \m^J 


Q^em > f 


(—) -/(O) 


where ovk = “em = e^/(47r) and the function f{z) is dehned as 


f {z) = 2 [ {1 + x)\og(^{l — x)z‘^ + x'^'^dx, 


for the quintuplet fermion and 


1 /* ^ 

/ {z) = —- / ^6 (1 — x) z^ + 9x‘^ — 4a; — 4j log^(l — x) z"^ + x'^'^dx, 


(5) 

( 6 ) 

(7) 


for the septuplet scalar. The septuplet scalar not only has SU{2)l gauge interactions, but 
also scalar couplings \^k and in the potential. These couplings also give a correction to 
the scalar masses. However, the mass corrections for all the components of the septuplet 
are exactly the same and no mass difference is generated in this way. We note that, 
although it may seem that our results differ from those in Ref. [9], we have checked 
explicitly that our expressions in Eqs. (5)-(7) are consistent with those in this reference. 

Inspection of Eqs. (5)-(7) reveals some relevant features of the mass degeneracy within 
the MDM multiplets. For light multiplets (m^ -C mw,'mz) the corrections are negligible 
for quintuplets but can reach a few GeV in the case of septuplets. Then, as the multiplet 
mass increases, the splitting in both cases approaches a common value. In fact, in the 
rriy. —)■ oo limit one hnds 


rriQ - mgi 


{Q^ 



g'2) X 166 MeV, 


( 8 ) 


in both MDM scenarios (quintuplet and septuplet). Therefore, although the behavior of 
the loop function f{z) is different for low values, it is the same for 3> rnwyiTT^z, 
leading to a universal splitting in case of heavy MDM multiplets. This universality can be 
used to estimate the resulting splitting at the two-loop level. The mass splitting between 
the singly charged and neutral components of a triplet fermion was calculated at the 
two-loop level in Ref. [31], Ending the value 164.4 MeV in the limit of an infinitely heavy 
triplet. Given the universality of this limit, we expect the same conclusion to hold in 
our two MDM scenarios. Finally, we also note that the previous expressions imply that 
the neutral component is the lightest of the components of the MDM multiplet, thus 
becoming a viable DM candidate. 


2.3 Neutrino mass matrix 

The Lagrangian relevant for the generation of neutrino masses is 

Cu = - + H.C., (9) 
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where = Viak^) is the usual neutrino Dirac mass term. Assuming {jnD)ia ^ 

the light neutrinos acquire Majorana masses via the canonical type-I seesaw mechanism, 

nty —rnk)m~^mD ■ ( 10 ) 

The light neutrino mass matrix in Eq. (10) is diagonalized by the Pontecorvo-Maki- 
Nakagawa-Sakata matrix Dpmns as 

/ mi 0 0 \ 

^PMNS ^PMNS = ( 0 m2 0 1 . (11) 

\ 0 0 m3 / 

It is common to use the Casas-Ibarra parametrization to write the neutrino Yukawa 
coupling y'' as [32] 

y — ’ (^- 2 ) 

where = diag(y'mj^), = diag(y'mi'), R is an orthogonal matrix {R^R = 

RR^ = 11) and the extra phase of the Yukawa coupling has been absorbed into the fields. 

As we will see below, non-thermal production of DM requires right-handed neutrinos 
with masses > 5 TeV. Therefore, the Yukawa couplings of the right-handed neutrinos 
participating in the generation of the light neutrino masses will be of the order of y'^ > 10“® 
in order to obtain an appropriate neutrino mass scale. 

3 Non-thermal production of Minimal Dark Matter 

The most recent measurement of the DM relic density by the Planck Collaboration is 
= 0.1186 ± 0.0020 at 68% CL [20]. In this section we will study the implications 
of this measurement for our setup, where the lightest neutral component of the MDM 
multiplet is the DM candidate due to the accidental Z 2 symmetry. 

A DM particle with SU{2)l gauge interactions and a mass larger than the W boson 
mass would mainly annihilate into a pair of gauge bosons, typically implying a strong 
reduction of the relic density. Thus, in the case of the MDM scenario with a large 
dimensional representation [9], the DM mass has to be above several TeV in order to 
reproduce the observed relic density. Moreover, the components of the MDM multiplet 
are always required to be nearly degenerate. This enhances the co-annihilation cross 
sections and also leads the Sommerfeld correction to the cross sections. Therefore this 
implies even heavier DM particles, with masses above 10 TeV [21,33]. Even if the DM 
mass is below the W boson mass, efficient co-annihilation with the charged particles in 
the MDM multiplet is still at work, and the DM relic density would be too low. 

In the absence of the heavy neutrinos, the above problem would be present in our 
models. However, the picture is slightly changed due to the non-thermal production 
of DM in the out-of-equilibrium decay of the right-handed neutrinos. The lightest right- 
handed neutrino is expected to be produced by a freeze-in mechanism in the early universe 
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Figure 1: Examples of solutions of the coupled Boltzmann equations. 


if its Yukawa couplings are small enough [34-38]. Then, once produced, it mainly decays 
into two body hnal states, Ni —)■ huai The decay width into two body hnal 

states is computed as 

(2/V^)ii”^7Vi 


Fwi — 


Stt 


(13) 


where the masses of the gauge and Higgs bosons have been neglected.^ There are also 
subdominant three body Ni decay processes into the components of the multiplet, such 
as Ni —)■ and Ni —)■ mediated by the gauge bosons. Then the charged 

particles decay into x°. These three body decay processes occur due to the mixing between 
left-handed and right-handed neutrinos. The branching ratio of these processes, including 
all the components of the MDM multiplet, will be denoted as Br,!. The number of DM 
particles produced per Ni decay is 2 Br^^, since a pair of DM particles is produced in each 
decay of Ni due to the conservation of the accidental Z 2 parity. In our analysis we will 
consider Br^ as a free parameter. Although the minimal models discussed in this paper 
predict a too low Br,^ value, we will comment below on how to increase this parameter 
with a minimal extension. 

The evolution of the number densities of the Ni and species in the early universe 
is given by the coupled Boltzmann equations 


dn 


Ni 


dt 

dnx> 

dt 


SHhn^ = 




gN^m%^m^oT f tun 


= -{(^eSv) 


nlo 


Ki 


eq 2 

n 0 


TTlyA) 




2BT^TN^nNi 


(14) 

(15) 


Here z = m-x/T with the temperature of the universe T, = 2 is the number of degrees 
of freedom of A^i, Ki{z) is the second modihed Bessel function, is the equilibrium num- 

^Note that in order to simplify the notation we have decided to denote the Ni decay width into two 
body final states as Fat, . However, this should not be confused with the 7Vi total decay width, which 
would also include three body final states. 
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Figure 2: Contours satisfying the observed DM relic density in the rn^o-T plane for 
several hxed Br^. The lightest right-handed neutrino mass is hxed to = 20 TeV for 
the upper panels, = 10 TeV for the middle ones and mjvi = 5 TeV for the bottom 
ones. The left panels are for the quintuplet fermion DM case while the right panels are 
for the septuplet scalar DM scenario. 


ber density of and (cTeffu) is the thermally averaged effective DM annihilation cross sec¬ 
tion, including co-annihilation processes with degenerate particles. The dominant annihi¬ 
lation and co-annihilation channels are and 

W^* —)■ jW^. In our analysis we approximately include the Sommerfeld effect for the 
effective annihilation cross section by using the results obtained in Ref. [21]. 

The DM relic density is obtained by solving the Boltzmann equations. For the nu- 
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merical analysis, it proves convenient to replace the number densities utvi and n^o by 
Yni = uni/s and Y^o = n^o/s, where s is the entropy density. Similarly, the time variable 
t is rewritten in terms of the temperature T. Some examples of numerical solutions are 
shown in Fig. 1, obtained for the quintuplet fermion DM scenario, with 2Br^ = 0.001, 
= 5 TeV and = 2 TeV. For each fixed DM mass, two different values of the 
decay width Fcan satisfy the observed relic density. For a small decay width F, the 
amount of produced Ni is moderate and they slowly decay after DM freeze-out (left panel 
in Fig. 1). For a large decay width F^v^, a lot of iVi’s are created and they decay somewhat 
fast into DM particles (right panel in Fig. 1). In this case, since the DM particles are still 
coupled with the thermal bath, the DM production due to the decay of Ni and the DM 
annihilation compete. The correct relic density is finally obtained. 

The parameter space satisfying the DM relic density measured by Planck for several 
hxed values of Br^ is shown in Fig. 2. One can see from the figure that for the quintuplet 
fermion scenario, the decay width of the lightest heavy neutrino Fshould satisfy 

IQ-i® GeV < F^, < 10-^2 GeV for tun, = 20 TeV, 

10"^® GeV < F^, < 10-^3 GeV for mjv, = 10 TeV, 

10"^® GeV < Ftvi < 10"^^ GeV for tun, = 5 TeV, 

in order to get the measured DM relic density with 2Br^ < 0.001. Thus, one finds 
that the required size of the lightest right-handed neutrino Yukawa coupling satisfies 
10”^^ ^ 2/^ ^ 10“®. For the septuplet scalar DM scenario, a larger branching ratio is 
required because the effective annihilation cross section is larger than in the quintuplet 
fermion DM case. In both cases we find that a much lower DM mass, compared to the 
standard MDM scenario, is possible. This is the main result of our paper. 

A threshold at m^o ^ 9.4 TeV is observed in the upper left panel of Fig. 2. When 
the DM mass is above 9.4 TeV, the effective annihilation cross section becomes too small, 
leading to > 0.12. Since the Ni decay can only increase the DM relic density, 

the m^o > 9.4 TeV region is excluded in our scenario. One should also note that a lower 
bound on the decay width Fcan be derived from Big Bang Nucleosynthesis (BBN). The 
conservative bound for the Ni lifetime < 0.1 s, which corresponds to F^vj > 10“^^ GeV, 
should be taken into account not to affect the predictions of BBN [39,40].® 

Since the Yukawa coupling for the non-thermal production of DM is required to be 
rather small, the lightest heavy neutrino Ni does not play any role in the generation of 
active neutrino masses. As a result of this, the active neutrino masses and mixings are 
explained by the other heavy neutrinos N 2 and Y 3 , with the Yukawa couplings of the 
order of y'^ ~ 10 “®. 

®If the Ni lifetime was longer than 10^ s, the DM particles produced by the decay of Ni would give 
an additional contribution to the number of effective neutrino species without affecting BBN. This is 
possible since the DM particles are relativistic due to their large kinetic energy [41]. However, this is not 
the case in our scenario, where the Ni lifetime is shorter. 
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We briefly comment on DM particles lighter than the W boson mass. For DM masses 
below mw, the DM relic density can be non-thermally produced with 2Br^ > 0.01 if the 
mass of the lightest right-handed neutrino is < 300 GeV. However, there are two 
reasons to disregard this scenario. First, the required branching ratio into DM is too 
large and calls for an extension of the model, as explained below. And second, and more 
importantly, this DM mass range is completely excluded by the constraints from collider 
experiments as discussed in Sec. 4.1. 

Getting the required Br^ value 

In the minimal models discussed in this paper one expects a small Br^ value for 
m^o ^ miY, mz- The main reason is intuitive: three body decays are suppressed compared 
to two body ones since the mediators of the three body decays, the and Z gauge 
bosons, are much lighter than Ni. In the septuplet scalar case, the Higgs boson can also 
mediate these decays via the scalar coupling However, the branching ratio Br^ cannot 
be made as large as required to account for the DM relic density even if the maximum 
value allowed by direct detection experiments is considered. 

This problem can be easily solved by adding a new heavy mediator for the three body 
decays. The simplest one would be a real scalar a, with couplings to a pair of y helds 
{Ua'^x^XX or 2/crO'y^y depending on the variant of the MDM scenario considered). The 
real scalar a would mix with the SM Higgs boson h and, as a consequence of this, the 
lightest right-handed neutrino Ni would decay into Ni —)■ xx^a- One should note that 
although the term aNfNj is also possible, this coupling should be small enough so that 
A^i is not in thermal equilibrium and can be produced by the freeze-in mechanism. We 
hnd that the required value of the h — a mixing angle is sin a ~ 0.1 and the mass of the 
new mediator {H, the heaviest mass eigenstate resulting from h and a) is rriH ^ for 
Ua = 0(0.1). The coupling is also constrained by DM direct detection experiments, and 
we have checked that the value Ha = 0(0.1) is consistent with the current LUX bound for 
TeV scale DM. We also have checked that the coupling y„ = 0(0.1) does not affect to the 
DM relic density and indirect detection, which will be discussed in the following section, 
since the electroweak interaction of the MDM multiplet is dominant. The required mixing 
angle is perfectly consistent with the current measurements of the Higgs properties by the 
LHC [42,43]. Finally, such a real scalar could be related to the breaking of an extra 
symmetry existing at higher energy scales such as, for example, U{1)l or U{1)b-Li where 
L and B are lepton and baryon numbers, respectively. 

4 Constraints 

In this section we review the most relevant constraints in our scenario. These come 
from collider searches as well as from dark matter indirect and direct detection experi¬ 
ments. 
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Figure 3: LEP limits on quintuplet fermion (left) and septuplet scalar (right) scenarios. 
In the derivation of these limits we have taken photon polar angles in the 46° < 9 < 135° 
range and photon energies > 6 GeV. 


4.1 Collider constraints in low dark matter mass scenarios 

Let us briefly comment on collider constraints in our scenario. These are of course 
only relevant for low dark matter masses and, as we will see, they strongly restrict this 
possibility. 

First, we consider LEP bounds. Collider constraints derived from mono-photon searches 
at LEP turn out to be very strong for low dark matter masses. Multi-charged particles 
plus a photon are produced by the process 

e^e“ ^ Z*/j* ■ 

Then, the multi-charged particles decay into 

X±± ^ ^ X°7r^7r± , 

producing soft pions due to the mass degeneracy among the components of the MDM 
multiplet X- Therefore, if these pions are not seen due to their low energies, the resulting 
signal at LEP is mono-photon plus missing energy e+e“ Using DELPHI data [44, 

45], the authors of Ref. [46] used this idea to set limits on the suppression scale A of 
four-fermion contact interactions of the type C >4 ~ ^XA^e. 

Here, we translate the lower limits of the suppression scale A into upper limits on 
the total production cross section for —)■ XA7) where x denotes x^; X^^ (a^^^)- 

Although the energy distribution of the mono-photon in our case is not exactly the same 
as that obtained with the contact interactions, this estimate should provide a rough limit 
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in our scenario. For our analysis we consider photons with polar angles 6 in the range 
45° < 9 < 135° and energies > 6 GeV. These kinematical cuts are required to be 
consistent with the detection capabilities of the High Density Projection Chamber of 
DELPHI. Our results are shown in Fig. 3, both for the quintuplet fermion and septuplet 
scalar. As one can see from the hgures, a slightly stronger constraint is obtained for the 
quintuplet fermion. DM masses m^o < 90 GeV for the quintuplet fermion and m^o < 
79 GeV for the septuplet scalar are excluded. We point out that very similar bounds have 
been obtained in the recent Ref. [30]. 

Regarding LHC bounds, these have been recently analyzed in [47]. For the quintuplet 
fermion case, the lower bound m^o > 267 GeV has been obtained with the LHC running at 
y/s = 8 TeV and 20.3 fb“^ of integrated luminosity in ATLAS and 19.5 fb“^ in CMS [47]. 
This bound will be improved up to m^o > 668 GeV if the LHC does not hnd a signal 
with ^/s = 14 TeV and 3 ab“^. 

4.2 Gamma-ray constraints 

The DM annihilation channels — )■ W^W~ and 77 induce indirect detec¬ 

tion signals of DM. At present time, the cross sections for these two annihilation processes 
are drastically affected by the non-perturbative Sommerfeld effect due to the low kinetic 
energy of the DM particles. This typically leads to relevant constraints from observa¬ 
tions of gamma-rays coming from dwarf spheroidal satellite galaxies or the galactic cen¬ 
ter [48,49]. Note that even for the W~^W~ annihilation channel, high energy gamma-rays 
are generated in the decay of the W boson. 

Let us elaborate on the so-called Sommerfeld effect. As we have already discussed, 
the components of the MDM multiplets are naturally degenerate, and the mass scale of 
the DM particles is much higher than that of the gauge boson masses. In this case, the 
usual perturbative calculation for annihilation cross sections is not valid because long- 
range Coulomb-like forces which imply the Yukawa forces with small mediator masses 
among the MDM components distort the plane wave function of the incoming DM two- 
body state. Hence, the annihilation cross sections must be calculated non-perturbatively 
by taking into account the Sommerfeld correction [50]. This was pointed out for the 
hrst time for a wino-like neutralino DM scenario in supersymmetric models [33,51-54]. 
A similar calculation was performed for MDM models in Ref. [21]. In our models, the 
calculation is basically the same, except for the mass range we focus on. Therefore, we 
must proceed to the re-evaluation of the Sommerfeld enhanced cross sections, in order to 
be able to compare to the current bounds from indirect detection experiments. 

Our numerical analysis follows Ref. [55] and includes only the s-wave component of 
the cross section. In order to obtain the correction factor for the amplitude induced by 
the Sommerfeld effect, we must solve a coupled Schrodinger equation in the presence of a 
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potential generated by long-range forces, 

1 (filJi 


m^o dr'^ 


, rriyov^ , 

+ 


( 16 ) 


where is the wave function of the two-body DM state and Vij is the potential matrix. 
The indices i,j run as i, j = 1 — 3 for the quintuplet fermion and as i, j = 1 — 4 for the 
septuplet scalar. The wave function ^|J and the potential matrix V are explicitly given by 

8A - 4A -2B 0 


{r\x^^X ) 


ifj = 


{r\x'^X ) 
(r|x°X°) 
for the quintuplet fermion and 
/ (r|x+++x—) \ 


1/ = I -2B 2A - A -3\/25 I , (17) 

0 - 3^25 0 


V' = 


V ( 


r\x^X^) 


1 / = 


/ 


/ 18A-9A -W 0 0 \ 

-W 8A-4A -55 0 

0 -55 2 A-A -6^25 

\ 0 0 -6^2 0 y 


(18) 


for the septuplet scalar. Here A = aem/r + aw cos^ /r, 5 = A = 

166 MeV and r is the distance between the two DM particles. The non-trivial factor 
\/2 appears in some matrix elements of the potential due to the different normalization 
between the neutral and charged states. 

We are interested in the annihilation of the neutral state, which is governed by physics 
at short distances and thus described by the wave function at the origin, '0(0)- Therefore, 
as discussed in Refs. [21,55,56], in order to evaluate the Sommerfeld correction factor, 
we must solve the Schrodinger equation with linear independent boundary conditions for 
irregular solutions at the origin. For the quintuplet fermion, these boundary conditions 
at the origin and at r ^ oo are explicitly given by 


^p{0) = 





d'0(oo) 

dr 


iki 0 0 

0 ik2 0 

0 0 i/cs 


V’(oo), (19) 


where ki = y — m^oVij(oo)/2. After solving the Schrodinger equation for each 

boundary condition at the origin, the Sommerfeld factor matrix Aij is given by 


Aij — 




( 1 ), 


re 


—ikjr 


( 20 ) 


where the superscript j implies j-th boundary condition. The above discussion is straight¬ 
forwardly extended to the septuplet scalar. 

Finally, the absorptive parts describing the DM annihilations for the W~^W~ and 77 
channels are given by 


2/2 5 3\/2 \ 2 / 16 4 0 \ 

Tww = ^\ 5 25/2 15/^2 , F,, = ^ 4 10, 

\ 3x/2 I 5 /V 2 9 / ""xo \ 0 0 0 / 


( 21 ) 
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for the quintuplet fermion and 


^ww — 


T^O?W 

Tn?o 


( 9 

24 

33 

18\/2 \ 


( 

36 

9 

0 \ 

24 

64 

88 

48^2 

y 27 ra 2 ^ 

36 

16 

4 

0 

33 

88 

121 

66\/2 

’ m\ 

9 

4 

1 

0 

y i 8\/2 

"7 

00 

66^2 

72 ) 


1 0 

0 

0 

0/ 


for the septuplet scalar. The DM annihilation cross sections can be computed from the 
absorptive parts and the Sommerfeld factor A as 


CVy/W — 2 {^ATy/Y/ (TV^^ = 2(/ir„3t)33, (23) 

for the quintuplet fermion and 

(jvww = 2 [AVwwA^) = 2 , (24) 

for the septuplet scalar. 

Our numerical results for the quintuplet fermion and septuplet scalar into W^W~ 
and 77 are shown in Fig. 4. The DM relative velocity is fixed to u = 10“^. Since the 
actual DM velocity has a distribution, the cross section should be averaged over the DM 
velocity in a more sophisticated analysis. For the W^W~ channel, the upper bound on 
gamma-ray observations coming from dwarf spheroidal galaxies (assuming NFW prohle) 
is also shown [48]. For the 77 channel, the blue points show the H.E.S.S. limit for 
monochromatic gamma-ray lines coming from the galactic centre, assuming in this case 
a Einasto profile [49]. As one can see from the figures, the cross sections largely exceed 
the bounds except in several dips. These dips can be interpreted in analogy with the 
Ramsauer-Townsend effect, the scattering of low energy electrons by atoms of a noble 
gas [57,58]. The Ramsauer-Townsend effect is caused by the electromagnetic interaction, 
and the positions of the dips strongly depend on the mass difference A [57,58]. Due to 
this effect, the incoming two DM particles pass through unaffected by the potential when 
the DM pair has a hxed energy, thus drastically decreasing the annihilation cross section. 
When the mass difference A is as large as (T(10) GeV, the dips tends to disappear since the 
electromagnetic transition between the DM and the charged states becomes unefficient. 
The resonant behaviour in Fig. 4 is caused by the bound state of the two DM particles 
with zero binding energy [52,53,57,58]. 

For the quintuplet fermion, the dips appear at DM masses around 2 TeV and 7.5 TeV 
for the W^W~ channel. The positions of the dips are very close to those in the 77 
channel. This is the only possibility to evade the gamma-ray constraints simultaneously 
for W^W~ and 77 final states. For the septuplet scalar, the same thing occurs at around 
5.5 TeV and 7 TeV. 

Note that the gamma-ray constraints depend on DM profiles. If a cored DM profile 
such as Burkert or Isothermal is assumed, these constraints are relaxed. Concerning this 
uncertainty, two orders of magnitude looser bounds are expected at most. Moreover since 
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Figure 4: Cross sections for W~^W and 77- The DM relative velocity 

is fixed to u = 10“^. 


the positions of the dips due to the Ramsauer-Townsend strongly depend on the mass 
difference A, if a larger mass difference, like 0{1) GeV, were possible by an extension of 
the model, a broader DM mass range would be able to satisfy the gamma-ray constraints. 

We point out that for the septuplet scalar DM, there are additional contributions to 
the annihilations into W~^W~ and 77 due to the scalar couplings and A^. These 
contributions would lead to an increase. However, they are found to be subdominant 
compared to the gauge boson loops (with the Sommerfeld correction calculated above). 

One should note that the above discussion of the dips due to the Ramsauer-Townsend 
effect is controversial. More detailed analyses have recently been done in Refs. [59,60]. 
According to these references, in addition to the W~^W~ and 77 channels, gamma-ray 
production in the ZZ, Z7 and 'W^W~x channels should also be taken into account. If 
these channels are included, the dips disappear and thermally produced MDM with a 
mass m^o 9.4 TeV for the quintuplet fermion or m^o ~ 25 TeV for the septuplet scalar 
is excluded in the case of cusp DM prohles such as NFW and Einasto, but still allowed 
for cored prohles like Isothermal.® For lower DM masses, as in our case, the gamma-ray 
constraint becomes stronger than in the thermal MDM scenario. There is, however, a 
valid region in the parameter space (with a specihc DM mass) when cored DM prohles 
are considered. 

4.3 Direct detection 

In the MDM scenario, elastic scattering with quarks is induced at the one-loop level 
via SU{2)l gauge interactions. Our setup with non-thermal DM production allows for 
lighter DM particles than in the in the standard MDM scenario with thermally produced 
DM. With non-thermal production we hnd valid DM masses below 7.5 TeV in both cases, 

® Another way to evade the gamma-ray constraint would be to assume a sub-dominant MDM scenario. 
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whereas in the standard thermal scenario one has 9.4 TeV for the quintuplet fermion 
and 25 TeV for the septuplet scalar. For this reason, we expect that the constraints 
derived from direct detection experiments will be stronger than those obtained in the 
usual thermal scenario. 

The one-loop spin independent elastic cross section with a proton was found in Ref. [9] 
to be about cXp ~ 10“^^ cm^. This value for the cross section may seem too large compared 
to the current experimental bound obtained by LUX for a DM mass m^o < 7.5 TeV. How¬ 
ever, according to recent calculations including two-loop diagrams including DM and glu¬ 
ons, the Higgs mass measured at the LHC and recent lattice simulations for the strangeness 
content of the nucleon, it turns out that the elastic cross section gets reduced due to par¬ 
tial cancellations, leading to Up ~ cm^ [61,62]. This is below the current LUX bound 
and testable by the future direct detection experiments such as XENONIT. Given that 
the DM particles in our scenario are lighter than those present in scenarios with thermally 
produced MDM, the coming direct detection experiments will also test our setup. 

5 Summary and Conclusions 

In this paper we have discussed an extension of the SM with three right-handed neu¬ 
trinos and a large SU{2)l multiplet. The SU{2 )l multiplet is either a quintuplet fermion 
or a septuplet scalar. Despite imposing no additional symmetry, the lightest neutral 
component of the multiplet can constitute the DM content of the universe because of 
an accidental symmetry, as in the conventional MDM scenario. Furthermore, neutrino 
masses are induced by the canonical Type-I seesaw mechanism. 

However, unlike the conventional MDM scenario with thermally produced DM, in our 
setup the DM particles are non-thermally produced by the decay of the heavy neutrinos. 
This allows to lower significantly the DM mass and still be compatible with the observed 
DM relic density. Instead of DM masses as large as 9.4 TeV or 25 TeV, the DM mass in 
our non-thermal scenario can be as light as a few TeV. 

Finally, we have considered several experimental constraints in our scenario. First, 
we have discussed the possibility of a DM mass below the W boson mass, which is ex¬ 
cluded due to strong constraints coming from mono-photon plus missing energy searches 
at LEP. In particular, we found that the mass ranges m^o < 90 GeV and m^o < 79 GeV 
are excluded for the quintuplet fermion and septuplet scalar cases, respectively. Next, 
we considered indirect detection constraints, especially relevant due to potentially large 
Sommerfeld enhancements. In fact, we found that the annihilation cross sections for the 
W^W~ and 77 channels are considerably affected by the Sommerfeld effect. The quin¬ 
tuplet fermion DM can evade the strong constraints of the gamma-ray experiments at 
around only m^o 2 TeV and 7.5 TeV due to the drastic decrease of the cross sections 
by the Ramsauer-Townsend effect. For the septuplet scalar DM, the same thing occurs 
and the DM is predicted to be 5.5 TeV or 7 TeV. In addition, if more detailed 
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analyses of the gamma-ray constraints in MDM scenarios are taken into acconnt, our 
scenario would be allowed only at specific DM masses for cored DM profiles. 
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A SU{2) Multiplet Notation 


There are two common ways to denote SU{2) multiplets: 

• Tensor notation: This is the usual choice, see Refs. [28,29]. A multiplet 4) would 

be represented by a symmetric tensor with some indices, (four indices for quintuplet 
^^jkl indices for septuplet where all indices can take values 1 or 2. 

• Vector notation: This is the choice made in this paper. In this case, the quintuplet 
and septuplet are simply represented by a vector of 5 and 7 elements. 


Since both notations are equally correct, this choice is just a matter of taste. In 
fact, a dictionary that translates the analytical expressions among notations can be easily 
found. Regarding the elements of the quintuplet and septuplet, the relation between the 
two notations is given by 





/ +$111111 \ 

+$1111 \ 


+ ^$111112 

_^2$1112 


+ yT5$111122 

-x/6<I)1122 

^ijklmn _ ^ 

-X/^$111222 

+2$1222 


_^$112222 

+ $2222 


+ \/6$122222 


^ _$222222 y 


(25) 


This allows us to write the quintuplet and septuplet y as shown in Eq. (1). 

We now comment on SU{2) septuplet direct products. The product 7 ® 7 can be 
decomposed as 


7 (g) 7 = 135 © 11 a © 95 © 7^ © 65 © 3^ © Is, 


(26) 
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where the indices S and A mean symmetric and anti-symmetric contractions. The sym¬ 
metry properties of septuplet contractions are fundamental in order to determine the 
number of relevant scalar couplings. For example, when one considers four kinds of 
singlets are obtained since anti-symmetric parts vanishT 

7(g) 7(g) 7® 7 

= (13 © 9 © 5 © 1) ® (13 © 9 © 5 © 1) 

D 1 © 1' © 1" © 1'". (27) 

However one can check that only two of them are linearly independent. 
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